Quantum holonomies with Josephson-junction devices 
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We examined properties of a Josephson-junction system composed of two coupled Cooper-pair 
boxes (charge qubits) as a candidate for observation of quantum holonomies. We construct a uni- 
versal set of transformations in a twofold degenerate ground state, and discuss the effects of noise 
in the system. 
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I. INTRODUCTION 

Appearance of the geometric phases in physical sys- 
tems has been, since their first systematic treatment by 
Berry ^ , among the most fascinating physical phenom- 
ena. Generalized to systems with degenerate spectrum, - 
they give the possibility to devise systems in which the 
dynamical contribution is only the overall phase factor, 
and the actual transformations are of purely geometric 
origin. 

Quantum geometric transformations (below referred to 
as holonomies) have attracted even more attention af- 
ter Pachos et al. '' proved their potential use in quan- 
tum computing. As compared to the ordinary, dynamical 
computation, quantum gates are here realized by cyclic 
evolution of parameters, and the result depends only on 
the geometry of the traversed path. 

Construction of holonomies has been discussed for vari- 
ous physical realizations, among them are also supercon- 
ducting nanocircuits."^' ' In these realizations the neces- 
sary number of independent tunable parameters results 
in high complexity of the considered systems. Here, mo- 
tivated also by the recent experiment with Josephson- 
junction system composed of two coupled charge qubits,'' 
we consider a similar design as a potential candidate for 
observation of quantum holonomies. As compared to the 
proposals in Refs."'' ', where the simplest two-dimensional 
holonomies are constructed using four coupled charge 
qubits, we achieve substantial simplification using only 
two qubits. Also, we realize the transformations within 
a twofold degenerate ground state, rather than excited 
state (as in^'^). In this way we avoid the problem of de- 
population of the subspace in which the holonomies are 
realized. 

We begin with describing the system and defining the 
operational subspace. In the four-dimensional space of 
two charge qubits we find a configuration of parameters 
for which the ground state is degenerate. We perform 
the holonomies by selecting two of the parameters and 
varying them adiabatically along a certain loop. During 
the transformation the degeneracy of the ground state is 
maintained by adjusting the remaining parameters as a 
function of the others. 

The scheme presented here may be realized in a system 
without strong constraints on its parameters provided 



that the noise level is low. However, since the system 
of charge qubits is usually affected by the charge fluctu- 
ations, in the following wc optimize the design in order 
to suppress the noise (for certain models of errors we 
find a decohercnce-frec subspace'''^'''). Finally we discuss 
possible extensions of the scheme. The specific system 
considered here is described by the Hamiltonian of two 
coupled qubits. In the discussion we also comment on the 
possible application of our scheme to an arbitrary system 
with this model Hamiltonian. 



II. THE SYSTEM 

The system we consider consists of two "charge 
qubits". They are coupled to each other via tunable 
Josephson junction - a symmetric dc-SQUID (supercon- 
ducting quantum interference device) (see Fig. 1). As 
discussed later, we need to have the possibility to switch 
off all the Josephson energies completely (and we use 
dc-SQUIDS instead of simple junctions). On the other 
hand, to construct nontrivial transformations, the ampli- 
tude and phase (or equivalently the real and imaginary 
parts) of the couplings Ji and J2 should be controllable 
during the operations. To achieve this we replace one of 
the junctions in the left and right SQUIDs with a further 
dc-SQUIDs (the same technique has been used in Ref.^). 

If the dimensionless gate charges '^gi(2) = 
Cgi(2)^i(2)/2e are close to 1/2, the only relevant 
charge states for each island are |0) and |1) (zero or one 
extra Cooper pair on the island). The Hamiltonian may 
be then reduced to the four-dimensional space and in 
the charge basis |ni,n2) has the form 
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Here the diagonal elements arc the electrostatic 
terms -B„i,„2 = Eci{ngi - ni)^ + Ec2{ng2 - ^2)^ -t- 
Em{ngi — ni){ng2 — 71.2) [ni{n2) is the number of 
excess Cooper pairs on the first (second) island] with 
Ed = 4e2Cs2/2(CsiCs2 - C^J, and similarly for E,2. 



Here C- 



Sl(2) 



is the sum of all capacitances connected 



to the first (second) island. The electrostatic coupling 
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FIG. 1: The system of two coupled charge qubits used to con- 
struct the simplest non-Abelian holonomies. The supercon- 
ducting islands ("Cooper-pair boxes") are denoted by black 
nodes, ni and n2 are the numbers of excess Cooper pairs on 
the islands. The Josephson coupling between the two islands 
may be turned off completely. The couplings Ji and J2 can 
be on one hand turned off completely, on the other hand we 
may obtain complex couplings by changing the fluxes $ci(2)- 



the operations they are kept constant). 

The operational subspace for the simplest non-AbeUan 
holonomies is spanned by two lowest-energy eigenstates. 
To make the transformations of purely geometric nature 
we need to make them degenerate. This property in our 
system does not follow from the symmetry arguments 
and we need to control it by making the parameters not 
fully independent. In the following we will compare the 
energies of the lowest states, and from this constraint 
calculate one of the parameters as the function of the 
remaining ones. In this way we define the degeneracy 
domain in the control manifold. 

For the following configuration of the parameters 



Jrn = Jl= J2 = 0, Ugi = ?lg2 = 1/2, 



(2) 



the Hamiltonian Eq.(l) is diagonal in the charge basis, 
and the ground state is twofold degenerate and spanned 
by the states 



term due to finite capacitance of the middle dc-SQUID 
equals E,n — 4e^Cm/(CsiCE2 — C,^„) (wc do not assume 
weak electrostatic coupling and the capacitance Cm does 
not need to be small). The Josephson coupling J^^^ = 

\JiJjl^ - Jfr^? + ^J^jl^jfr C0s2(^<I.i) CXp[-*7/;($i)] 

(and similarly for J*-^'), where tanV'(^'i) = 
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jj^Van7r$i.^^ Here J. 



the Josephson energy of the junctions in the dc-SQUIDs. 
The term Jj> = Jjr{a) costt^c is the tunable Josephson 

couphng of the small SQUID. Finally J,„ — J™^ cos7r$m 
is the Josephson energy of the middle SQUID. Here the 
fiuxes are in units of $0 = hc/2e, the superconducting 
fiux quantum. We assume that the system is operated 
in the charge regime {Ec^2),Em » Jj^(r)' -^Mr)^ ■^^'')- 

A similar system has been used already in the context 
of geometric phases. Falci et al. used two coupled 
charge qubits to devise an experiment in which the Berry 
phase could be detected. Also, as already mentioned, the 
experiment by Pashkin et al. '' shows that a system of 
two coupled charge qubits can be well controlled and ma- 
nipulated with available experimental techniques. The 
difference between our system and those used in Refs.'"'^-^ 
is that we couple the islands using dc-SQUID instead of 
a capacitor, and we require controllable Josephson cou- 
plings Ji(2). 

The set of controllable parameters defines the control 
manifold. The Hamiltonian is characterized by the tun- 
able terms: Ed, Ec2, Ji, J2, Jm- We will use in the 
discussion rather the parameters that can be controlled 
directly, i.e., {^i,^2,^7n,ngi,ng2). The fluxes $ci(2) are 
used only to control the symmetry of the SQUIDs, and 
we will not consider them to be used as the actual tun- 
able parameters in the following discussion (between the 
operations they are kept at the values for which the cou- 
plings Ji(2) real, and may be tuned to zero. Just 
before the operations we can switch them to different 
values for which the dc-SQUIDs arc asymmetric. During 



We will refer to this particular point in the parameter 
space as the starting point, or the "no-op" regime (as 
this configuration will be kept between the sequences of 
transformations), and to the states |0) and |1) as the 
logical basis. 

Before we perform any transformations the system is 
initialized to the state |0) by turning off all the Josephson 
couplings, and tuning the electrostatic energy Eqi to be 
the lowest one. After long enough time the system will 
relax, and we switch to the starting point. We perform a 
holonomic transformation by varying the parameters adi- 
abatically along a suitable loop. During the operations 
the other charge states may be involved as well. The 
system, however, at any instant of time remains in the 
twofold degenerate ground state, and the resulting uni- 
tary transformation (holonomy) is limited to the logical 
basis only. The correspondence between the loop, and 
the performed unitary transformation is given by 



Ur = V cxp 



(4) 



where P is the path-ordering operator, F the path, R = 
{Xi} is the set of parameters (control manifold), and the 
2x2 matrix Ai is the Wilczek-Zee connection,^ 



Af = {a{R)\^\(3{R))- 



(5) 



I a) and |/3) are parameter-dependent states spanning the 
operational subspace, and for the 'no-op' regime they can 
be either |0) or |i). 

To construct the holonomies we need to find a loop 
corresponding to the desired transformation. Any uni- 
tary transformation can be factorizcd as a product of 
two noncommuting rotations in the qubit space. Here 
we find two such holonomies, namely, exp(i(f>iaz) and 
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exp{i(j)2crx). For our purposes it is convenient to interpret 
these rotations as phase shifts between |0) and |1), and 
between |+) = (|0) + |1)) /^/2 and |-) = (|0) - |1)) /V^, 
respectively. For given rotation, at each point of the con- 
structed paths the Wilczek-Zee connection A is then di- 
agonal [in the basis |0(1)) and | + (— )) respectively], and 
we omit the path-ordering operation. 



III. CONSTRUCTING THE HOLONOMIES 

The holonomies we find are in our case the aforemen- 
tioned rotations, e*'^i'^= and e*'^^'^==. The first transfor- 
mation wc consider is the phase shift e^'f'^'^' . If we keep 
Ji = Jm = (and change the flux <i>c2 to make the cou- 
pling J2 complex), the Hamiltonian simplifies to 



Hi = 
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(6) 



The obtained block form may be solved for each part in- 
dependently. We calculate the desired phase shift as the 
difference between the Berry phases for each block. At 
the same time, to make the phase shift of purely geomet- 
ric origin, we compare the energies of the lower states 
within each block to find the degeneracy subspace in 
the control manifold. Since the energies of the lowest 
states are functions of ngi,ng2, and $2, from the condi- 
tion Ea{ngi,ng2,^2) - Ef3{ngi,ng2,^2) = we find the 
degeneracy condition in the functional form ngi{ng2, $2)- 
This dependence is shown in the upper part of Fig. 2 

The parameter-dependent spectrum (for the degener- 
acy domain in the control manifold) is shown in the mid- 
dle plot of Fig. 2. The effectively four-level system has 
only three distinct energy values due to the degeneracy 
of the ground state. 

Finally the phase shift may be presented conveniently 
as the integral, 

01 = -^ / dng2d<t>2[Va{ng2,'i>2)-Vi3{ng2,<^>2)], 
^ Js{r) 

_ (7) 

where ^(r) is the domain in the ng2-^2 plane limited by 
the path F, and the "Berry field" is given by' 



q(/3) 



{aiP)\VRHi\ea{eb)) x (ea(efc)|Vflgi|a(/3)) 

{Ea{/3) — Eea{eb))'^ 



(8) 

Here jca) and |eb) are the excited states from the upper- 
left and the lower-right block of Hi respectively, and 
Eea{eb) thcir energies. The field difference V = Va — Vp 
is shown in the lower plot of Fig. 2. 

Similarly, we can find the second generic holonomy, 
namely, e^'^^"'^ . For Ji = J2 = J (the coupling should 
be complex during the operations, so the fluxes ^ci{2) 
should have different values than at the starting point, 
but chosen in this way that the Josephson energies may 
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FIG. 2: Construction of the exp{i(l>iaz) rotation. Upper plot 
- the degeneracy subspace calculated from the comparison of 
the energies of the lowest-energy levels; middle plot - the 
parameter-dependent spectrum of the system (the ground 
state is twofold degenerate); lower plot - the "Berry field". 
All plots taken for the following values of system param- 
eters: E^i/E^ = 0.8, £c2/-E„ = 1, J^/E,^ = 0.1, and 
Jf/E,„ = 0.05. 



be equal), Eqi — Eiq = E, the Hamiltonian may be 
again presented in the block form, this time in the basis 
{|00),|ll),|±) = (|01)±|10))/^/2}: 



Ho = 
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Comparison of the lowest energy of the upper-left block 
with E + Jml^ gives the degeneracy subspace, as shown 
in the upper plot of Fig. 3. The middle plot shows the 
parameter-dependent spectrum, again the ground state 
being twofold degenerate. Since the state |— ) does not 
vary with the parameters, the desired phase is simply the 
Berry phase acquired by the ground state in the upper- 
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left block. The phase here equals 

h^l [ dng2d^V'ing2,<i>). (10) 
^ Js{r) 

The field V' is shown in the lower plot of Fig. 3. 



r 




FIG. 3: Construction of the exp{i(j>20-x) rotation. Upper plot 
- the degeneracy subspace calculated from the comparison of 
the energies of the lowest-energy levels; middle plot - the 
parameter-dependent spectrum of the system (the ground 
state is twofold degenerate); lower plot - the "Berry field". 
All plots taken for the following values of system parame- 
ters: Eci/E„, = 0.8, Ec2/E^ = 1, j'^^V-Em = 0.1, and 
Ji''^'' / Em ~ 0.05. Here we assumed for simplicity that the dc- 
SQUIDs are identical so that Ji — J2 implies $1 = $2 = ^■ 
In practice only the first condition needs to be satisfied 



IV. DISSIPATIVE EFFECTS 

The scheme described above is valid for a system with- 
out any strong constraints on its parameters - the only 
idealization so far is the assumption of perfect symmetry 
of the small dc-SQUIDs. As already noted, the opera- 
tions are performed within a ground state, and the re- 
laxation does not depopulate the operational subspace. 



However, as the degeneracy is maintained by external 
parameters, we should expect that the system will suf- 
fer from their fluctuations. Now we will give qualitative 
analysis of the noise in the considered system, and de- 
scribe how the design can be optimized in order to sup- 
press the dissipation. In Joscphson charge qubits the 
Ohmic noise and the low-frequency charge fluctuations 
coming from the substrate arc usually relatively strong. 
Since the system is operated in the charge regime, and 
both sources couple to the charge on the islands we apply 
here the spin-boson model with the system Hamiltonian 
approximated by its electrostatic part only 

iOi 1 uji 2 T 1 2 t 

k 

+ {al + al) X + alYi + a^,Y2. (11) 

Here X, Yi are the bath operators. We divided the inter- 
action term into the part describing the common envi- 
ronment for the two islands, and part with independent 
environments. If now the two islands have equal energy 
splitting {uji = UJ2 = ^o)^ a-nd the dominant interaction 
with environment can be modeled by common environ- 
ment, the subspace {|01), |10)} is well protected (in the 
limiting case of completely correlated noise = the 
states span the decoherence-free subspace). Our goal is 
then to minimize the uncorrelated noise in the islands as 
much as possible. To achieve this for the noise coming 
from the substrate we can, as discussed in Ref. ' ' place 
the islands close to each other. To suppress the Ohmic 
noise we could use a common voltage source controlling 
the bias of the islands (as in the experiment by Pashkin 
et al. ''). As we present here, this may work well also for 
our purposes, provided that the islands are nearly iden- 
tical. We use again two voltage sources, but now their 
roles are not symmetric. Instead, the main, common 
source (the fluctuations of which give correlated noise) 
will be used to tune the system, while the second, auxil- 
iary, is connected to one island via very small capacitance 
(compared to Cgi(2)) as it is used to tune a small (as ex- 
plained below) difference between the gate charges Ugi 
and 71^2 • The resulting noncorrelated noise will be then 
much lower than the correlated noise of the main source. 

To be more specific, the fluctuations in Ug should be 
correlated (rather than in Vg) and we require that the 
electrostatic properties of both islands are nearly identi- 
cal {]Eci -Ec2\= SEc < £^01(2), and Cgi - Cg2 = 5Cg < 
Cgi(2))- Then the common voltage source Vgi = Vg2 
gives first of all nearly equal energy splitting of the is- 
lands [lui « UJ2 in Eq.(ll)], and also relates the param- 
eters Ugi w ng2. We note that for the starting point as 
well as during the operation e"^^'^^ we have £'01 = £10. 
This for identical islands gives ngi = ng2- Also during 
the operation e*'^^'^- there is a region in the parameter 
space in which the relation Ugi = ngi{ng2,^2) nearly 
coincides with Ugi = ng2 [this region is determined by 
Eqi = Eio < Eoo,Eii, where the states |0) and |T) have 
equal and lowest electrostatic energies. The small dif- 
ference between ng2 and ngi{ng2,^2) is in this region 
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a result of small, but nonzero Josephson energies]. We 
then can use an additional voltage source that is coupled 
only to the first island. This voltage will give an addi- 
tional term Sugi which compensates the effect of small 
differences in the properties of the islands, and during the 
gj0io-= operation the small difference ng2 — ngi{ng2,^2)- 
Since each gate can be performed using various paths, 
we may even further improve the fidelity of the operations 
by examining the sensitivity of the degeneracy splitting 
to the fluctuations of each parameter, and choose paths 
from the region least affected by the fluctuations. For 
instance, close to the starting point [defined in Eq.(2)] 
the main voltage Vg fluctuations to the first order do not 
shift the system out of the degeneracy subspace. The 
same applies to $2 for the e"^'^"'"' operation and to <& 
for e*'^^'^^. This can be easily seen in the upper plots 
of Figs. 2 and 3. Those fiuctuations lie within a plane 
tangential to the degeneracy subspace. To protect the 
system also during the operation we should construct the 
loops relatively close to the starting point. 

V. DISCUSSION 

To summarize, we have found that the system proper- 
ties give us the opportunity to realize any unitary trans- 
formation within a qubit space using gates of purely ge- 
ometric origin. Assuming the perfect performance not 
affected by the noise we may rely on the simplest de- 
sign without strong constraints on the system parame- 
ters. However, such a system composed of two solid- 
state qubits is usually strongly probed by the environ- 
ment which results in fast dephasing. Our qualitative 
analysis made for this particular system specifies the con- 
ditions under which the system should be well protected 
even in the presence of noise. 

We have not discussed here the coupling between such 
elementary blocks, as the question of applicability of 
holonomies in quantum computing remains open and re- 
quires detailed analysis of resulting fault tolerance for 



individual proposals. However, as the logical basis here 
is defined by states differing by the charge configuration 
of the islands, such extension based on the method de- 
scribed in Ref.'' should be adjustable to our system. 

The scheme presented here may be applied to any 
system described by the Hamiltonian of two interacting 
qubits of the form 

H = Bi-a^ +B2-(7^ + JzCfUl + J_i(cr|cr2 + h.c), (12) 

provided the parameters are tunable. Since there are un- 
til now many well-developed schemes concerning qubit 
control as well as the controlled qubit-qubit coupling 
(usually for the sake of quantum computation), realiza- 
tion of quantum holonomies may appear very natural 
and straightforward in such systems. For the same rea- 
son the model of a four-dimensional quantum system as 
a candidate to construct and observe holonomies may 
turn out to be more applicable than the simplest, and 
widely used in implementation-independent discussions 
three-level model (see for example '). 

As far as the Josephson-junction systems are consid- 
ered, we may also use the device to study the process 
of the adiabatic charge transport: ' the logical basis is 
spanned by the states with one excess Cooper-pair on 
one island and the second with zero. Quite naturally the 
second discussed operation, generated by Ux, corresponds 
for (p2 = 7r/2 to the charge pumping cycle. ' 
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